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Abstract 

Nucleon-antinucleon annihilation into a near backward- (or forward-) produced meson and a 
high invariant mass lepton pair admits a factorized description in terms of antinucleon (or nu- 
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cleon) distribution amplitudes (DAs) and nucleon to meson (or antinucleon to meson) transition 
distribution amplitudes (TDAs). We estimate the cross section of backward (and forward) pion 



and r/-meson production in association with a high invariant mass lepton pair for the kinemat- 



ical conditions of GSI-FAIR. The cross sections are found to be large enough to be measured 
with the PANDA detector. Interesting phenomenological applications of the approach are thus 



expected. 



PACS numbers: 13.60.-r, 13.60.Le, 14.20.Dh 



1. INTRODUCTION 

The ambitious experimental program for the PANDA facility at GSI-FAIR (see jl[ 
will provide the experimental access to new classes of hard reactions for which the 
factorized description applies. Study of these reactions will open a new eyehole to look on 
hadron's interior and will bring a new perspective to the problem of hadronic structure 
in terms of the fundamental degrees of freedom of QCD. 

A tempting possibility is to consider nucleon-antinucleon annihilation into a lepton pair 
with a high invariant mass q 2 = Q 2 in association with a light meson Ai = {n, rj, p, u, ...}: 

N(Pn)N( Pn ) -> 7*(q)M(p M ) -+ i + (pe + )r( Pe -)M( PM ). (1) 

The analysis of jjj, Q], identifies two similar factorization regimes for the reaction ([1]) 
corresponding to forward and backward peaks of the meson production cross section. 

Let us now specify our conventions 1 . The z-axis is chosen along the colliding nucleon- 
antinucleon. For the PANDA setup it is natural to select as the positive direction of z 
the direction the antinucleon is moving in the NN center-of-mass (CMS) frame. 

• Within these conventions, the kinematical regime in which the w-channel momentum 
transfer squared, u = (pn ~Pm) 2 , is small corresponds to the meson moving in the 
direction of the initial nucleon that is the backward direction. 

• Analogously, the kinematical regime in which the t-channel momentum transfer 
squared, t = (p^ — P m) 2 , is small, corresponds to the meson moving in the direction 
of the initial antinucleon that is the forward direction. 

A detailed account of the corresponding kinematics is presented in Appendix |A] 

The backward and forward factorization regimes for the reaction (CQ) are casted with 
the help of two collinear factorization theorems schematically depicted on Fig. [TJ 

1. For backward kinematics the collinear factorization theorem (see the left panel of 
Fig. [T]) is valid once s = (pn +Pn) 2 = W 2 and the invariant mass of the lepton pair 
q 2 = Q 2 are large; the corresponding skewness variable £ M ( 1A2j) . which characterizes 



1 Definitions of t and u here match Ref . [5[ and differ from those of Rcf . 
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the w-channel longitudinal momentum transfer, is fixed; \u\ is small as compared to 
Q 2 and W 2 . The meson Ai is produced in the near backward direction in the NN 
CMS. We refer to this kinematics as the backward, or the w-channel, factorization 
regime. 

2. Analogously, the t-channel collinear factorization theorem presented on the right 
panel of Fig. [1] is valid once W 2 and Q 2 are large; the corresponding skewness 
variable ( lAllj) is fixed; \t\ is small as compared to Q 2 and W 2 . The meson Ai is 
produced in the near forward direction in NN CMS, and we refer to this kinematics 
as the forward, or the t-channel, factorization regime. 

It is worth emphasizing that for the moment the rigorous proof of the collinear factor- 
ization theorems for meson production in association with a high invariant mass lepton 
pair in NN annihilation as well as for backward meson electroproduction does not exist. 
Some physical arguments in favor of this kind of factorization for the spacelike version of 
the reaction ((Jj) were presented in Refs. 6|, lZJ. In this paper we do not aim to prove the 
relevant factorization theorems but rather perform the feasibility estimates to try to see 
some evidences of the factorized description in the experiment. 

The factorization theorems for the reaction (TjQ) presented in Fig. [1] involve two kinds of 
nonperturbative objects: the conventional antibaryon (or baryon) distribution amplitudes 
(DAs) and the baryon to meson (or antibaryon to meson) transition distribution ampli- 
tudes (TDAs). The baryon to meson (antibaryon to meson) TDAs are defined through 
baryon (antibaryon)-meson matrix elements of the nonlocal three quark (antiquark) op- 



erator on the light cone 
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0^(X 1 n, A 2 n, X 3 n) = e clC2C3 ^ a (X in )^(X 2 n)^(X 3 ny, 
O a /3 7 pT X (Ain, A 2 n, X 3 n) = e ClC2C3 ^^ p (Ain)^ 2 T (A 2 n)^ 3 x (A 3 n). (2) 

Here a, (3, 7 stand for quark (antiquark) flavor indices and p, r, \ denote the Dirac spinor 
indices. Antisymmetrization stands over the color group indices Ci.2,3 2 . The gauge links 
in (j2J) can be omitted in the light-like gauge A + = 0. For the case when the meson Ai 
is a pseudoscalar, to the leading twist-3 accuracy, the form factor decomposition of AiN 



2 See App. IB II for the details of our index conventions. 
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FIG. 1: Two possibilities for collinear factorization of the annihilation process 
NN — > 7*(q , )A^(px). Left panel: backward kinematics (\u\ ~ 0) . Right panel: 
forward kinematics (\t\ ~ 0). N(N) DA stands for the distribution amplitude of 
antinucleon (nucleon); MN(MN) TDA stands for the transition distribution amplitude 
from a nucleon (antinucleon) to a meson; CF and CF' denote hard subprocess 

amplitudes (coefficient functions). 

(MN) matrix elements of (J2D involves eight invariant MN (MN) TDAs H MN (H mR ), 
each being a function of the three longitudinal momentum fractions Xi (Yli x i = of 
the corresponding skewness variable £ and of the momentum transfer squared A 2 , 



H( MN ^ MM \x 1 ,x 2 ,x 3 ,e'\A u ' t2 ) 

_ { y(MN), (MN) ^ A (MN), (MN)^ T (MNJ, (MN )}{ ^ ^ ^ ^ ^ 



(3) 



The isotopic symmetry and the charge conjugation invariance further reduce the number 
of independent nucleon to meson TDAs. In particular, C-invariance allows us to express 
all MN TDAs through MN TDAs (see App. EJ). 

Extensive studies of the properties and of the physical interpretation of ttN TDAs 



are presented in Refs. 



13Ml5l|. Conceptually, the A^A^ TDAs share many common 



features both with the GPDs and the nucleon DAs. Indeed, the crossing transformation 
relates A4A^ TDAs with the MN generalized distribution amplitudes (GDAs), defined as 
the matrix element of the same light cone operator between A^A^ state and the vacuum. 
In the pion case, this allows to establish useful constraints from the chiral dynamics since 
it N GDAs reduce to combinations of the usual nucleon DAs in the soft pion limit. In 



addition, irN TDAs were recently estimated within the light-cone quark model |16|]. The 
hadronic matrix elements of QCD operators are extensively studied on the lattice. In 
particular, irN matrix elements of the operators corresponding to the Mellin moments of 
ttN TDAs have been recently considered on the lattice by several groups, see [l7], 
and Refs. therein. 



On the other hand, similarly to the GPD case [1914221] . a comprehensible physical 
picture may be obtained by Fourier transforming M.N TDAs to the impact parameter 
space. 

2. MODELLING irN TDAS 

First estimates of the cross section of tt° production in association with a high-Q 2 
lepton pair in pp annihilation for PANDA@GSI-FAIR conditions within the factorized 
description involving ttN TDAs were presented in [4|. In this analysis, an oversimplified 
model for nN TDAs based on the soft pion theorem was employed. This model has 
obvious drawbacks: its domain of applicability is limited to the immediate vicinity of the 
soft-pion threshold; TDAs in this model do not have any intrinsic A 2 - dependence and 
finally, the cross channel leading baryon exchange contribution, which turns out to be 
dominant in many cases, was not taken into account. 



In Ref. [5], a two component model for ttN TDAs was proposed. It includes the spectral 
part casted in terms of quadruple distributions and a D-term like contribution which is 
determined by the nucleon pole exchange in the cross channel (see Fig. [2]). Quadruple 
distributions are fixed with the help of the soft pion theorem for tcN TDAs in terms of 
nucleon DAs. Once the pion mass is neglected (m = 0) in the strict soft pion limit (£ = 1, 
A 2 = M 2 ), this model gives the same predictions for the NN — > £ + £~tt cross section as 
the model employed in |4j. 
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FIG. 2: Nucleoli pole exchange model for ttN (left) and irN (right) TDAs; dashed circles 
contain typical LO graph for the nucleon (antinucleon) electromagnetic form factor in 
perturbative QCD; rectangles contain the nucleon pole contribution into ttN (ttN) TDAs. 

As it was found out in the contribution of the spectral part dominates the cross 
section for large values of £ ~ 1 (which is difficult to reach in experiments), while for 
small and intermediate values of £, the nucleon pole exchange part provides the dominant 
contribution to the NN — >■ 7*7r cross section. 

In Sec. 14 1[ we present our estimates of NN — > £ + £~tt cross section within our factor- 
ization scheme in the w-channel factorization regime 3 (near backward pion production) 
using the nucleon pole exchange model for ttN TDAs. 

The explicit expression for the contribution of the cross-channel nucleon exchange into 
the isospin-i ttN TDAs was obtained in 15|, 
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where V p , AP and T p stand for the nucleon DAs; g 7 



NN 



13 is the pion-nucleon phe- 



3 In the remainder of the paper, in most cases we omit superscripts u and t for the kinematical quantities 
(£, A 2 , p, n etc.) referring to the factorization regime. 
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nomenological coupling and 

3 

e E RBLOl, X 2 , X 3 ) = JJ 6{0 < X k < 20 (5) 

fc=l 

ensures the pure Efremov-Radyushkin-Brodsky-Lepage (ERBL) support. 
For 7i~n TDAs, we get 

(ir~ n) 



JV(940) 



{V 1>2 , A 1>2 , T± 2} (xi, £, A 2 ) = V2{V li2 , A 1>2 , T li2 }^(x 4 ,£,A 

L ' JV(940) k J 



{T 3 ,T 4 }^ nJ (* 4 ,£,A 2 ) =0. (6) 

k J 7V(940) 



We also present some estimates of t 
component model for irN TDAs of Ref. 



re NN — > £ + £ 7T cross section within the two 

I 



3. CROSS SECTION OF NN -> 7*vr -> ir£+£- 

Below we review the formulas for the cross section of 

N{p N , s N ) + N{ PN , s R ) -> 7 *(g) + vr( P7r ) -> + r (p/-) + tc{ P7T ) (7) 



established in [4j]. The starting point is the general formula for the unpolarized differential 



cross section of the reaction (see e.g. [231 ]) 

1 



" g= 2(2,)»A(^,M 2 ,AP) '- M ^- | ^< LIPS) - (8 > 

where A(x,y, z) = \J x 2 + y 2 + z 2 — 2xy — 2xz — 2yz is the usual Mandelstam function. 

Following the standard procedure, the 3-particle differential Lorentz invariant phase 
space rf 3 (LIPS) is decomposed into two 2-particle phase subspaces: those of 7*7r and £ + £~ 
systems. The former can be easily computed in NN CMS while the latter is computed 
in £+£- CMS yielding the result, 

da ~ 2(2tt)5A(^ 2 ,M 2 ,M 2 ) | - M 1 8W 2 [ ^ 8 ' 1 j 

where = d cos #* dtp*^ in NN CMS. By rffi^ = dcos9ed<pe, we denote the produced 
lepton solid angle in £ + £~ CMS. By expressing cos6** through u = (j>n — P-k) 2 , 

du = ^^A(W 2 , M\ M 2 )A(W 2 , Q 2 , m 2 ) (10) 



and integrating over the azimuthal angle 92* of the produced pion and over the azimuthal 
angle of the lepton pair (p#, the following formula for the differential cross section of the 
reaction (jTJ) is established: 



- J Y 1 1 11 



dudQ 2 d cos 9 e 64W 2 (W 2 - 4M 2 )(2vr) 4 ' 



The next step is to express the average-squared amplitude \J\/l NN ^ i+ ^^\ 2 through the 
amplitude A4 s ) ^ Sn of NN — > 7*7r. Within the factorized approach of [4|, at leading order 
in a s , the amplitude J\4 s ^ Sn reads 

1 



S s x » s *x(^ A z ) - S'™ T{£, A 2 ) , (12) 



where 



n _ , (47ra fl ) 2 V47ra em /^ 

Here, f n = 93 MeV is the pion weak decay constant and fx is a constant, which determines 
the value of the dimensional nucleon wave function at the origin 121 ]; a em ~ is the 



electromagnetic fine structure constant and a s ~ | is the strong coupling constant. 
The spin structures S and S 1 are defined as 

s s»s N = v(p R ,s N )t(X) l5 U( PN ,s N ); 

= Lv{p Si SR )t(\)A Tl5 U(p N , s N ). (14) 

We employ the standard "hat" notation for the contraction of a four-vector with the Dirac 
matrices. e(A) stands for the polarization vector of the virtual photon. I and V denote 
the convolution integrals of itN TDAs and antinucleon DAs with the hard scattering 
kernels computed from the set of 21 relevant scattering diagrams (see [3]): 

{1, !'}(£, A 2 ) = / dx l dx 2 dx 3 8(x l + x 2 + x 3 -2£) 
J-i+z J-i+z J-i+i 



: fdyxf dy 2 [ dy 3 5( Vl + y 2 + y 3 - 1) (2 V{i? a , R'J + R'j] . 

J o Jo Jo \ a=1 a=8 ) 



(15) 



The integrals in x^s (y^s) stand over the support of ixN TDA (antinucleon DA). 

For the w-channel factorization regime of pp — > ir°£ + £~, the coefficients R a and R' a 
(a = 1, .., 14) correspond to the coefficients T a and T' a presented in Table I of [3] up 
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to an overall irrelevant phase factor ^jv -1 ?^ 3 and to the replacement — is — > is in the 
denominators of hard scattering kernels. This latter change mirrors the fact that we 
consider 7* in the final state rather than in the initial state as it is for backward pion 
electroproduction. 



We also note that Ref. 
TDAs, which does not maintain the po 
paper, we use the parametrization of 



13] employed a somewhat inadequate parametrization of ttN 
ynomiality property in its simple form. In this 
lj which maintains the property of polynomial^. 
The relation between the two parametrizations is presented in App. A of [5|. For np — > 
ir~£ + £~ reaction the coefficients R a and R' a are the same as for pp — » tc°£ + £~ channel up 
to an obvious change 

{V 1>2 , A 1>2 , T li2i 3, 4 } ( " 0p) -> W, 2 , A 1>2 , T li2i 3, 4 } { "" n) . (16) 

The differential cross section for unpolarized beam and unpolarized target is computed 
from the squared amplitude ( T12]) averaged over spins of the initial particles: 

\M^\ 2 = \j2 M T s HMrn* (ir) 



with the given combination of photon helicities A, A' = T, L. \M.ll\ \M-lt\ and |A^tl| 
vanish at the leading twist accuracy. We employ the relation 

XyCA^A) = -gT + -^—(j/W+pW) (18) 



(p ■ n) 

A7 1 



to sum over the transverse polarizations of the virtual photon and we get 

1 2 2(1 + Q 
4 1 1 £Q 6 11 1 M 2 



\m t \^y, ra 2 = i^i 2 ^^(m 2 - UpT), (19) 



A j 1 

where C is defined in (fi3|) . 

The averaged- squared amplitude for the process (I7j) reads 

\M»S^*-«\ a = \ £ ^rie 2 Tr{^e(A)p £+ e*(A')}^(.Mr P T- (20) 

As pointed above, at the leading twist accuracy, only the transverse polarization of the 
virtual photon is contributing. After computing the relevant trace in ( 120]) in the £ + £~ 
CMS and integrating over the lepton polar angle (p#, we finally get: 

|2 27re 2 (l + cos 2 9 t ) 



d<p e \M NN ^ e+£ -*\ 2 = \M T \ 2 v ' 2 • (21) 

Leading twist LJ 
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A kinematical cut in Ay, that is equivalent to a cut in cos 8* (1A7|) . is a convenient 
way to select near backward and near forward kinematical configurations. On Fig. [3] with 
solid lines, we show the dependence of the CMS scattering angles 8* (see eqs. (1A7I) , 
( 1A16I) ) for the w-channel (left panel) and the t-channel (right panel) factorization regimes 
as the function of A^ in and A^ in respectively, for W 2 = 5 GeV 2 , £ = 0.3. For these 
kinematical parameters, the A^ in = —0.2 GeV 2 cut corresponds to selecting pions flying 
backward into a cone with the aperture of ~ 120° in the NN CMS. 

In order to implement the effect of the Al| cut, we consider the cross section integrated 
over the appropriate bin in u. Indeed, u can be expressed as follows as the function of £ 
and A 2 -: 

_ a 2 .^ + 1) 2e(M 2 (i-p-m 2 (e + i)) , 99 , 
u ~ i-c i-e ' 1 ] 

Therefore, the bin Ay min < A 2 - < corresponds to the following bin in it: u min < u < 



Mmax, where 



„ {P A 2 , A^ + l) , 2^(M 2 (l-Q-m 2 (e + l)) . 
2e(M 2 (l-Q-m 2 (e + l)) 

Wmax(0 = " ( 23 ) 
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u-channel regime; W 2 =5 GeV 2 ; £"=0.3 t-channel regime; W z =5 GeV 2 ; f'=0.3 
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FIG. 3: The solid lines illustrate the dependence of the CMS scattering angles 6* for 
the w-channel (left panel) and the t-channel (right panel) factorization regimes for the 
reaction ([7]) as functions of A^ in and A^ 2 ^, respectively. The dashed lines illustrate the 
dependence of the LAB frame scattering angles 9^ AB for the two factorization regimes as 
the function of A^ min and A Tmin , respectively. Note that the forward peak is narrowed 
and the backward peak is broadened due to the effect of the boost from the CMS to the 
LAB frame which corresponds to the nucleon N at rest in the PANDA setup. 



4. PSEUDOSCALAR MESON PRODUCTION IN ASSOCIATION WITH A 
HIGH INVARIANT MASS LEPTON PAIR IN NN ANNIHILATION 

1. Results for 7r° and 7r production 

On Figs. H] and IU we show our model predictions for the integrated cross section 
da f" m » r da 

dQ^^L iu r i ^m^e t (24) 

of pp — > £ + £~tt° and of of pn — > £ + £~tc~ as the function of Q 2 , for several values of W 2 
(W 2 = 5, 10 and 20 GeV 2 ) with the cut at A^ min = -0.2 GeV 2 . As phenomenological 
input, our model for ttN TDAs requires the nucleon DAs V p , A p , T p at the normalization 
scale fi 2 ~ Q 2 . ^A vas t literature exists on the phenomenological solutions for nucleon 



DAs (see e.g. 24j, |25] for the discussion). In order to quantify the sensitivity of our 
model prediction on the input nucleon DAs, we show the cross section estimates for 
the case of several phenomenological solutions fitting the nucleon electromagnetic form 
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factor: Chernyak-Ogloblin-Zhitnitsky (COZ) 26] 4 (dashed line with long dashes), King 
and Sachrajda (KS) |27j 5 (solid line), Braun-Lenz-Wittmann (BLW NLO) model of [25| 6 
(dashed line with medium dashes) and NNLO modification of BLW model suggested 



in Ref. 28| 7 (dashed line with short dashes). The phenomenological solutions for the 
nucleon DA presented in the literature are usually given at a low normalization point 
fi 2 = 1 GeV 2 , while the adequate choice of the normalization scale within our factorized 
approach is fi 2 = Q 2 . In our cross section estimates presented on Figs. HI [5] we do not 
take systematically into account the evolution effects which are expected to be not very 
significant due to the small lever arm in Q 2 . 

To get an estimate of these effects for the kinematics in question we show on Fig. [6] 
the ratio of the integrated cross sections 

da/dQ 2 \^ =Q 2 



da/dQ 



21 



(25) 



\/i 2 =l GeV 2 

for pp — > tt°£ + £~ with and without account of the evolution effects for input nucleon DA 
as the function of Q 2 (with a s kept fixed). The result is presented for the BLW NNLO 
input nucleon DA. For Q 2 = 2 4- 4 GeV 2 the effect of the evolution results in a mere 40% 
decrease of the cross section. This effect is much below the model dependence due to 
different phenomenological input. For other input phenomenological DAs we expect the 
effect of the evolution to be of a similar order of magnitude. 

Finally we conclude that within the considered parameter range, the cross sections of 
pp — > tt°£ + £~ and pn — > tt~£ + £~ presented on Figs. HJ E] seem to be large enough for a 
detailed investigation to be carried at PANDA@GSI-FAIR. 



4 Eq. (13) of Ref. [2a. 

5 Eq. (4.6) of Ref. [231 

6 See Appendix B of Rf 

7 See Table I of Ref. m 
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pp -> n f + r- W 2 =5 GeV 2 ; \A 2 T \< 0.2 GeV 2 ; 



pp -> TrVr; W 2 =10 GeV 2 ; |A||< 0.2 GeV 2 





2.5 3 3.5 4 2 2.5 3 3.5 

Q 2 [GeV 2 ] g 2 [GeV 2 ] 

pp -> A + r; W 2 =20 GeV 2 ; |A^|< 0.2 GeV 2 ; 




2.5 3 3.5 

Q 2 [GeV 2 ] 



FIG. 4: Integrated cross section da/dQ 2 for pp — )■ £ + £~7r° as a function of Q 2 for different 
values of W 2 = 5, 10 and 20 GeV 2 for various phenomenological nucleon DA solutions: 
COZ (long dashes); KS (solid line); BLW NLO (medium dashes) and NNLO modifica- 
tion [28] of BLW (short dashes). 



It is worth mentioning that the cross section that we obtained with the BLW NLO 
model input is much smaller than that in the case of models including next-to-next-to- 
leading conformal spin contribution (such as COZ, KS and NNLO modification of BLW). 
In fact, the pp — > 7r°7* and np — > 7r _ 7* amplitudes turn to be zero at the leading twist 
accuracy once we employ the asymptotic form of the nucleon DA as input. This zero is 
reminiscent of the zero for the proton electromagnetic form factor with the asymptotic 
DA (lOj. The cross section thus turns out to be small for the BLW NLO as well as for 



the Bolz-Kroll (BK) 29( input nucleon DAs which are known to be rather close to the 
asymptotic form of the nucleon DA. 



13 



pn -> n~tC; W 2 =5 GeV 2 ; \A 2 T \< 0.2 GeV 2 ; 



pn -> 7T/ + r; W 2 =10 GeV 2 ; |A||s 0.2 GeV 2 
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p« -> iCtT; W 2 =20 GeV 2 ; |Af |< 0.2 GeV 2 
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FIG. 5: Integrated cross section da/dQ 2 for — )■ £ + £~7r~ as a function of Q 2 , for 
different values of W 2 = 5, 10 and 20 GeV 2 for various phenomenological nucleon DA 
solutions: COZ (long dashes); KS (solid line); BLW NLO (medium dashes) and NNLO 
modification [28| of BLW (short dashes). 



pp -> TiVr; BLWNNLO input 




3.0 
Q 2 [GeV 21 



FIG. 6: The ratio of the integrated cross sections ^ 2 q2 

account of the evolution effects function of Q . 



with and without taking 
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In Fig. [7J we show the differential cross section da /dQ 2 d cos 8* for pp — > £ + £~tt° as a 
function of cosd*, for W 2 = 5 GeV 2 and Q 2 = 2.5 GeV 2 , both for the near forward and 
backward factorization regimes. As a consequence of C-invariance, the cross section of 
([T]) within the t-channel factorization mechanism can be obtained from that within the 
it-channel factorization mechanism with the change ( 1C1I) (see discussion in App. IB 31) . 
The COZ solution for the nucleon DAs is used here as numerical input. The forward 
and backward peaks produced within the suggested factorization mechanism are clearly 
visible. They look perfectly symmetric in the CMS frame. However, it is worth mentioning 
that when boosting from the CMS to the LAB frame (which corresponds to the nucleon 
N at rest in the PANDA setup) , the forward peak is narrowed and the backward peak is 
broadened (see Fig. [3]) . 

pp -> tt° ( + r ; W 2 =5 GeV 2 ; Q 2 =2.5 GeV 2 ; 

^1000 
> 

3. 800 
&h 600 



b 

T3 200 



-1.0 -0.5 0.0 0.5 1.0 

cos 0* 

FIG. 7: Differential cross section da/dQ 2 dcos9l for pp — > £ + £~tt° as a function of cos6** 
for W 2 = 5GeV 2 and Q 2 = 2.5GeV 2 . Forward and backward peaks are clearly visible. 
COZ solution for the nucleon DAs is used as the numerical input. Dotted region denotes 
scattering over large angles in which the present factorization description does not apply. 

Obviously, the ttN TDAs in the nucleon pole model have purely ERBL support and 
turn into zero at the borders of the ERBL domain. In a general model in which 7rN TDAs 
do not vanish at the borders of the ERBL domain (so-called crossover trajectories w = — £, 
v = ±£'), important contributions into the cross section may come from the convolutions 
with coefficient functions that are highly singular at the crossover trajectories. 

To estimate this effect we consider the two component model for 7riV TDAs [5] which 
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includes the spectral part (formulated in terms of quadruple distributions) and the nucleon 
pole part, which is the analogue of the D-term familiar from the GPD case. The spectral 
part is nonzero in both the ERBL and the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP)-like domains and does not vanish at the cross over trajectories. This model is 
an attempt to maximally take into account the general requirements (following from the 
underlaying fundamental field theory), such as the following: 

• polynomiality property and the support property which are the consequences of the 
Lorentz invariance; 

• permutation symmetry and the isotopic invariance; 

• chiral properties manifest through the soft pion theorem. 

Nevertheless the two component model of Ref. 5| is still very flexible. In particular, the 
factorized Ansatz for the quadruple distributions and the shape of the profile function 
are ad hoc assumptions. So we can rather make a model dependent estimate of the 
contributions of the cross over trajectories and of the DGLAP region into the cross section. 
As the nucleon pole contribution dies out for £ = 1 ,the spectral part contribution becomes 
relatively more important for higher £. 

In Fig. [8] we show the effect of adding the contribution of the spectral part for the 
integrated cross section da/dQ 2 of pp — > tt°£ + £~ presented inn Fig. HI Thin dashes show 
the pure nucleon pole part contribution; thick dashed lines show the sum (at the amplitude 
level) of two contributions. The COZ phenomenological solution for the nucleon DA is 
used as the numerical input. For different phenomenological inputs the effect of adding 
the spectral part is similar. 
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FIG. 8: The integrated cross section da/dQ 2 of pp — > £ + £~tc° as a function of Q 2 for 
VT 2 = 5 GeV 2 (left panel) and W 2 = 10 GeV 2 (right panel) in the two component model 
of ttN TDAs (thick dashed line) v.s. the nucleon pole model (thin dashed line). COZ 
nucleon DA is used as the numerical input. 



One may see that for W 2 = 5 GeV 2 , adding the contribution of the spectral part results, 
at most, in a factor of 2 increase of the cross section for Q 2 ~ 4 GeV 2 (this kinematical 
region corresponds to rather large values of £ ~ 0.6). This effect is not very important 
in view of the large uncertainty due to different input phenomenological DAs. Therefore 
the account of the spectral part contribution together with the question of appropriate 
modelling of quadruple distributions will become a topical question after a stable signal 
from the backward (and forward) regimes will be available from PANDA. Moreover, for 
W 2 = 10 GeV 2 the Q 2 interval in question corresponds to rather small values of £ ~ 0.1, 
for which the nucleon pole part totally dominates over the spectral part within the present 
model. 



2. Results for r\ production 

In an analogous way, we may consider near backward (and forward) r] meson production 
in association with a lepton pair 

N(p p ) + N(p p ) -> 7 *(g) + v(Pv) ~+ + t~(Pt-) + v(Pv)- (26) 

The nucleon pole model for r]N TDAs is similar to that for ttN TDAs (jlj) with the 
obvious change of phenomenological coupling g^NN ~ * g^NN- Estimates of g n NN and 
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dr/NN phenomenological couplings taken from Table 9.1 of Ref. (301 give ^ ~ 0.3. The 

9-irNN 

formulas for the pp — > 7*77 — > £ + £~r] are obtained from the relevant formulas of Sec. [3] 
with the obvious change of values of the masses and the couplings. In Fig. [9j we show 
the results for the integrated cross section da mt /dQ 2 for pp — > £ + £~7] as a function of 
Q 2 , for different W 2 = 5, 10 and 20 GeV 2 for various phenomenological nucleon DA 
solutions: COZ (long dashes), KS (solid line), BLW NLO (medium dashes) and NNLO 



modification 



of BLW model (short dashes). 



pp -> 7] tC; W 2 =5 GeV 2 ; |A||< 0.2 GeV 2 ; 
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FIG. 9: Integrated cross section da/dQ 2 for pp — > £ + £~i] as a function of Q 2 for different 
values of W 2 = 5, 10 and 20 GeV 2 for various phenomenological nucleon DA solutions: 
COZ (long dashes); KS (solid line); BLW NLO (medium dashes) and NNLO modifica- 



tion 



of BLW (short dashes). 
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5. CONCLUSIONS 



In this work, we have estimated the cross section for pseudoscalar meson production 
in association with a high invariant mass lepton pair in nucleon antinucleon annihilation. 
We have performed this analysis in the framework of the QCD collinear factorization. 
We used the model for the TDAs which was introduced recently in [ijj] and for which 
the major input is the parametrization of the nucleon distribution amplitudes. For most 
of the choices, the magnitudes of the cross sections are sufficient to be measurable with 
the luminosity foreseen at the PANDA experiment at GSI-FAIR. On the other hand, it 
turns out that the cross section depends much on the parametrization used for these DAs, 
which may be analyzed in a positive way as an indication that the experimental study 
of these processes will help to distinguish between the different models for these DAs. 
Note that these different models all claim to be in reasonable agreement with the nucleon 
form factor measurements. We are, however, still far from proposing a strategy for a 
model-independent extraction of the TDAs from experimental data. More observables, 
and in particular spin asymmetries, should be discussed. Coefficient functions should also 
be calculated at next-to-leading order and the effect of the Q 2 evolution of the TDAs 
should also be studied in details. Note that these NLO order effects are known to be 



different in timelike and spacelike processes 



3l|. Meanwhile, the preparation for the ex- 



perimental analysis requires detailed background studies (such as that of the final states 
where a 7r + 7r~ meson pair mimics the lepton pair) and a careful estimation of detector 
efficiency. Although we concentrated here on pseudoscalar (tt and 77) meson production, 
this framework can be extended to the production of vector mesons (p, u) which will 
uncover complementary aspects of hadronic structure both experimentally and theoreti- 
cally. Also a purely hadronic description of pp — > is proposed in Refs. [32| and 
331 ]. respectively, in a Born model and a Regge model description of baryonic exchanges. 
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A. KINEMATICS 

Below we review the relevant kinematical quantities discussed in J4] for the two possible 
factorization regimes for the reaction (CQ) involving nucleon to meson and antinucleon to 
meson TDAs. 



n-channel factorization regime 

Within the conventions natural for the PANDA setup, the u-channel factorization 
regime (see the left panel of Fig. [1]) in which \u\ is small corresponds to the meson moving 
in the direction of the initial proton (backward direction). 

The z-axis is chosen along the colliding nucleon-antinucleon. We define the light-cone 
vectors p u and n u such that 2p u ■ n u = 1. We introduce the average momentum and 
u-channel momentum transfer 



P u = \{Pm+Pn)\ A u =p M -p N ; A 



u2 



U. 



(Al) 



A u t is the transverse component of A u (A^ 2 < 0). We define the skewness variable 

A" • n u 



2pu . n u 

The following Sudakov decomposition is valid for the momenta of the reaction (jj]) 

M 2 „ 



(A2) 



Pn 
Pn 



C)p u 



-n 



2M 2 (1 + C) 



W 2 - 2M 2 + WVW 2 - AM 2 

2 _ \u2 

pm = (i-Op u + ' 



W 2 _ 2M 2 + WVW 2 - AM 2 B 
p H — n u - 



2(1 + e) 



A" = -2Cp u + 



M 2 



m 



- AVp 2 



-n 



M 2 



i + e 



2c + ^(i+n + o(i/w 4 ) 



n u + Aj.; 

W 2 — M 2 



p u + 



m 



A^ 



T 



1 - 



+ 0(1/W 2 



n u - A^, 



(A3) 
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where M (m) denotes nucleoli (meson) mass. 

In the nucleon rest frame (which corresponds to the PANDA laboratory frame) the 
light-cone vectors p u and n u read 

M 1 + f u 

P U \N rest = 2(I+e) {1 ' 0, Q '~ 1}; ^-t = ^^{1,0,0,1}. ( A4 ) 

With the help of the appropriate boost we establish the expressions for the light-cone 
vectors in the NN CMS 

P u \nn cms = K, 0, 0, -a u }- n u \ RN CMS = {(3 U , 0, 0, ?>}, (A5) 

where 

a - 4(1 + e) ; ^ ~ • (Abj 

The meson A4 scattering angle in the NN CMS for the w-channel factorization regime 
then can be expressed as: 

-(i - eK + V^/ 3 " 

cos^ = ; — = - — ■ (A7) 



One may check that for A^ 2 = indeed cos 6^ = —1, which means backward scattering. 
We also quote some useful relations for the kinematical quantities: 

M 2 m 2 



i + e i - e 



1f u AM 2 
Q 2 = Q 2 = f^W 2 + u-3M 2 + ^ + 0(1/^ 2 ); 

1 — f u AM 2 
t=(PM~ Pn) 2 = ~Y^ W2 + ™ 2 ~M 2 + ^ + 0(1/^ 2 ). (A* 



This allows to express £ M as 

5 2W 2 -Q 2 + u-3M r y ' 

i-channel factorization regime 

The t-channel factorization regime (see the right panel of Fig. [T]) in which \t\ is small 
corresponds to the meson moving in the direction of the initial antiproton (forward direc- 
tion). 
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We define the light-cone vectors p l and n* such that 2p* • n* = I. 

P* = \{pm + Pn); ^ = Pm-Pn\ A t2 = t. 

t 2 



(A10) 



AV is the transverse component of A* (A^, • A^ = Atp < 0). We define the skewness 



variable 



A* • n* 



2P t ■ n* ' 

The following Sudakov decomposition is valid for the momenta of particles: 

M 2 t 



PN = (l+C)p t + 



-n 



Pn 



2M 2 (1 +(*) 



W 2 - 2M 2 + WVW 2 - AM 2 



W 2 _ 2M 2 + W\/W 2 - AM 2 , 

p + , 7* n ; 



2(1 +e) 



PM = (l-e)p t + m2 1 _f 2 n t + A t T ; 



+ 



M 2 



m 2 - Atp 2 M 2 



i-e i+e 



n* + A^; 



2^ + ^(l+a + 0(l/^ 4 ) 



V + 



W 2 — M 2 m 2 - A* 2 



(All) 



i+e 



+ 0{l/W 2 ) 



n l - A* T . 



(A12) 



In the antinucleon rest frame the light-cone vectors p 1 and n l are 



P I Wrest — 



-{1,0,0,1}; n'U rcst = — f {1,0,0,-1}. 



(A13) 



2(1 + e) L "'"'"' ~ J ' " IJvrusL 2M 
The explicit expressions for the light-cone vectors p* and n* in the NN CMS read: 

p'Inncms = {«*, 0,0, a*}; n% N CMS = {£*, 0, 0, -£*}, (A14) 



where a and /r are defined as 



t _ W + VW 2 -4M 2 _ t _ - V ^ 2 - AM 2 ) (1 + £*) 



(A15) 



4(1 +£*) ' ^ 4M 2 

The meson A4 scattering angle in the iVA CMS for the t-channel factorization regime 
then expresses as: 



2 -A' 2 



i - ew - n -^hp 



cos 6 



M 



((i-ey-^?-^ 2 



(A16) 
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Again one may check that for A^, 2 = 0, indeed, cos 6^, = 1, which means forward scat- 
tering. 

We also work out the relations: 

M 2 m 2 



K 2 = ^4 ( t - 2£ 



i + e i-e. 

IP 1 4 M 2 

Q 2 = q 2 = ^TT t W 2 + t - 3M 2 + i^- + 0(1/W 2 ); 
l + 4 1 + 4 

1 — A* 4/lf 2 
« = (PM ~ Pn? = -T-T t W 2 + m 2 -M 2 + —— + 0(1/W 2 ). (A17) 

1 + 4 { 1 + 4 

and express £*: 

t ^ Q 2 -t-M 2 
4 ~ 2W 2 -Q 2 + t-3M 2 ' 

B. CHARGE CONJUGATION AND G-PARITY ISSUES 



(A18) 



1. Notations and conventions 



Below we employ our usual system notations of Ref. 

• Letters from the beginning of the Greek alphabet are reserved for the SU(2) isospin 
indices a, (3, 7, 1, k = 1, 2. 

• We have to distinguish between upper (contravariant) and lower (covariant) SU(2) 
isospin indices. We introduce the totally antisymmetric tensor e a /3 for lowering 
indices and e al3 for rising indices (E12 = s 12 = 1). 

• Letters from the middle of the Greek alphabet A, fi, v denote the Lorentz indices. 

• Letters from the second half of the Greek alphabet p, r, x are reserved for the Dirac 
indices. 

• Letters from the beginning of the Latin alphabet a, b, c ... are reserved for indices of 
the adjoint representation of the SU(2) isospin group. 

• Letters ci, C2, C3 stand for SU(3) color indices. 
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The nucleon field N a (N a ) transforms according to the covariant (contravariant) rep- 
resentation of the isospin SU(2). We adopt the following standard convention for the 
nucleon field [34]: 

J^T E {e tkx d^(k,s)V(k,s) + e- tkx b a (k,s)U(k,s)}; 
j2^fY E {e ikx tf a (k,s)U(k,s) + e- ikx d a {k,s)V(k,s)}. (Bl) 

Here spinors U(k, s) and U (k, s) = U'(k, s)7o describe a nucleon, respectively, in the initial 
and final states, while spinors V(k, s) = V'(k, 3)70 and V(k, s) describe an antinucleon in 
the initial and final states. 

The creation and annihilation operators in (IBli satisfy the usual anticommutation 
relations for fermions 

{P(p, 8 ), b^(p',s')} = (2n) s ^5 3 (p-p>)6 ss ,5%; 

{d a (p,s), d^(p',s>)} = (27r) 3 ^(p-p')6 ss ,6P. (B2) 

The "in" nucleon state \N a ) is defined according to: 

|Aq) = \N p (p,s))=b{(p,s)\0)- \N 2 ) = \N n (p,s)) = bl(p,s)\0). (B3) 

Analogously, the "in" antiparticle state \N a ) is defined as: 

IN 1 ) = \N^p,s)) = d^(p,s)\0); \N 2 ) = \N n (p, a)) = £ 2 {p, a)|0> . (B4) 

The charge conjugation operator C has the following effect on the nucleon creation and 
antinucleon annihilation operators 



CV(p,\)& = VN d L (p,X); (B5) 

Cd.'faXp = VNbt(p,\), (B6) 

where r] N denotes the nucleon field charge parity. These equalities are to be understood 

not as equalities of SU(2) tensors but rather as equalities for their components. From 
(IB 6 1) and the hermitian conjugation of (1B5I) . we get, respectively: 

C\N") = VN \N,}; 

C\N L )= V * N \W). (B7) 
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2. Antinucleon DA 



The isotopic parametrization for the nucleoli DA reads [1 



4(01^^3^^(1)^(2)^37(3)1^^)) 

= e^Af^ 13 >(l, 2, 3) + e^f M^ 12 >(1, 2, 3), (B8) 

where the invariant isospin amplitudes read 

M£f >(1, 2, 3) = « x ^ p (l, 2, 3) + < x ^(l, 2, 3) + ^2^(1, 2, 3)) ; 
M£ X 13 >(1,2,3) = M p ^ T 12 >(l,3,2). (B9) 

Here {u^, a^, ^Jp^x are the conventional Dirac structures, 

< x = (PC)Mu{p)) x ; < x = (P7 5 ^)p.(f/W)x ; = (° P »c) P A^u(p)) x , 

(BIO) 



where C is the charge conjugation matrix; a^ v = |[7 M , 7"] and cr pi/ = p M cr 



[IV 



The effect of the charge conjugation operator on the quark field is 34| : 



C^ a tf = r) q CW a , (Bll) 

where transposition refers to the Dirac index. C stands for the charge conjugation matrix 
and 7] q is the corresponding charge parity. 

Now using (IBlip together with (IB7[) . we establish the link between nucleon and antin- 
ucleon DAs, 

4(0|e ClC2C3 ^-(l)^(2)^37 (3 )| i v tb ) ) 

= 4^(0|£ clC2C3 C t C^^"(l)C t CvI/^(2)C t CvI/^(3)C t |iV t (p)) 

= 4^ 3 (0|e ClC2C3 {C^) p {C% c >) r {C^) x |A»>. (B12) 

Now we can define the antinucleon DA as the SU(2) tensor: 

4(0k ClC2C3 {^f) r {*«) x \N\ P )) = e a ^M*[ 13 \l,2,3) + e a ^M^ 2 \l,2,3), 

(B13) 
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where 

M^ 12 >(1, 2, 3) = f N (v* x V*(l, 2, 3) + < x ^(l, 2, 3) + tf T ^(l, 2, 3)) ; 
i< x 13 >(l, 2, 3) = M%W(1, 3, 2), (B14) 

where the relevant Dirac structures are 

(t£ x ) r = (Ct) w ,(Ct) TT ,(Ct) xx ^ r , x , = (Cpf^X; 

(< X ) T = (c%Ac%Ac% x ^ T , y = {c lh p) T pr {v) T x 

(t% >x ) T = {^) pp ,{C ] )rAC\ x d N p , Tl . x , = {Ca^iV^l (B15) 
Lifting the transposition with respect to the Dirac indices in (1B15I) one gets 

< x = (Cp)p.(^7 5 ) x ; 

t% >x = (Ca m ) pT (VYls) x . (B16) 

Therefore, we conclude that as the consequence of charge conjugation invariance an- 
tiproton DAs occurring in (1B14I) are expressed as follows through the usual proton DAs 
of (|B9|): 

{V, A, TY( yi ,y 2 ,y 3 ) = -^{V, A, Ty( yi ,y 2 ,y 3 ). (B17) 

V N Vg 

3. Antinucleon to pion TDAs 

Below we establish the relations between the antinucleon to pion and nucleon to pion 
TDAs. Our main line of reasoning is analogous to that in the case of the antinucleon 
DAs. However, charged pions are not eigenstates of the charge conjugation operator C. 
To work out the isotopic formalism for the antinucleon to pion TDAs it is convenient to 
employ the concept of G-parity (see e.g. [35]), which is the combination of the charge 
conjugation and the rotation around the second axis of the isospin space, 

g = Ce inh . (B18) 

First, let us work out the effect of Q on the quark and antiquark fields 

Q^ a Q^ = Ce i7Th ^ a e- i7Th Cl (B19) 
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According to the conventions of [15(, we choose to transform ty a field according to the 
covariant representation and to transform field according to the contravariant repre- 
sentation of the isospin SU(2): 

[/., * J = I {<r a f a ; [/., * a ] = -\ (a a f a V a , (B20) 

where I a stand for the group generators and a a are the Pauli matrices. 
The effect of the isospin rotation around the second axis then reads: 

^(tfV*' 3 = ( cos + sin f H^) O^V (B21) 
Thus, we conclude that 

Gufl = v q (-l)(Cd) T p] QdpQ^ = r) q (Cu) T p . (B22) 

We also need to specify the effect of Q on pion "out" states and nucleon "in" states. 
The triplet of pions has negative G parity, 

(7r|0t = —(tt]. (B23) 

Employing our isospin conventions for the nucleon states: summarized in Appendix IB II 
we conclude that: 

Q\N P ) = V* N (-m n }; Q\N n ) = V * N \NP). (B24) 
For antinucleon to pion TDAs we can write down the following parametrization 
4(7r a |0^ 7pTX (l, 2, 3)|A>> 

= (/ fl )W M ^x )3/2 (l,2,3) + ^K)' 7 M^ )1/2{13} (1>2,3) 
+e aj (a a YpM^ )l/2 2, 3) , (B25) 

where the tensor {f a )\ a p 7 }i which is totally symmetric in a, (3, 7, reads: 

since (crj)i = (a a )J = M\ ^ (B26) 

The isospin and permutation symmetry identities for the isospin- 1 and isospin- 1 in- 
variant nN TDAs M (7rAr)l /2 i 12 ), M {nN)l / 2 {13} and M (7rAr)s / 2 are the same as those for the 
relevant ttN TDAs (see Sec. 4 of Ref. In particular, 

M^v )l/2{13} (1,2,3) = M (^)v 2 {m (M)2) _ (B27) 
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Now, using G-parity operator Q instead of C, we repeat the derivation of App. IB 2l and 
establish the link between ttN and ttN TDAs. 
For example, 

\/2 

-— M (wAr)3 /2(i ; 2,3) + >/2M (irJV) va(l,2,3) 
3 



7T 



- | ^ t ^ p (l)6;t^ T (2)^ x (3)^t^| i v n ) 



X 



(-l)<^(C) pp ,(C) TT ,(C) xx ,(7r-|(d P J (l)KO' (2)(%') i (3)|iV p ) 

XX 

{ - X ( M Xv /2 ^' 2 ' 3 )) T + ^ (M^ 2 (l, 2, 3)) T }. (B28) 



One may check that the general parametrization for nN TDAs (IB25j) is consistent with 



that for 7i N TDAs (see eq. (54) of Ref. [13]) once 



M^ 2 (l,2,3)) r = ^(Cl) P AC\AC% x ,M { ; T N ^ /2 '(1,2,3); 



M^ 2 (1,2,3)) T = ^^^^(C^^K^^xx'^ 37 ^!^^), (B29) 



where the transposition refers to the Dirac indices. The Dirac structures s nN = 
{V12, a i2> ^1^34} occurring in the parametrization of the isospin-| and isospin-| in- 
variant amplitudes M^'^ 2 and M^ nN ^ 3 ^ 2 are defined by 



> t 

N 



S pT,X 



) = {&) P( ,{tf) TT ,{CX x >sf T , )X „ (B30) 



where sffi, , are the Dirac structures occurring in the parametrization of irN TDAs (see 



eq. (12) of Ref. [13]). For the relevant nN TDAs we, thus, get 



H^( Xl ,x 2 ,x 3 ,^A 2 ) = -^—H^ N (x 1 ,x 2 ,x 3 ,^A 2 ). (B31) 

V N Vq 

C. AMPLITUDE CALCULATION WITHIN THE it-CHANNEL FACTORIZA- 
TION MECHANISM 

Below we argue that the amplitude of (CQ) within the t-channel factorization mechanism 
can be obtained from that within the w-channel factorization mechanism (see Sec. [2]) with 
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the obvious change of the kinematical variables: 



Pn -»■ Pn, Pn -> Pjv! 
A M — >■ A* («-►*); 

r -> (ci) 

Since isotopic invariance plays no particular role in the present consideration we consider 
one flavor of quarks and one "pion" with negative C-parity: 77^ = —1. Generalization for 
the physical case is straightforward with the use of the G-parity formalism described in 
App.E! 

We introduce the following parametrization for the corresponding matrix elements: 

{irlVpMVrMVxixJlNtpN)) = A pT (p N ,A u )(BU( PN )) x H(x t ,C,u) 
{0\^(yiM(V2)%{V3)\N(px)) = D- pf (p N )(V( PN )E)^(y t ), (C2) 

where A(p^,A u ), B, D(p N ) and E are the relevant Dirac structures; H(xi,£, u ,u) and 
$>(yi) denote the scalar TDAs and DAs respectively. 

We now express the matrix elements occurring in the t-channel factorization scheme 
through the Dirac structures of the u-channel factorization scheme occurring in flC2j) : 

(n\^ T p (x 1 )^(x 2 )^ T x (x 3 )\N(p ff )) 
'In'Iq 'h 

= -±j-^AC^) pT ^BU) x H{x l} e, t); (C3) 
'In'Iq 'h 



*M*f(y2)*M\ N M) 

1 rtC&C^C&(0\^ yi )%(y 2 )^y 2 )\N(p N )) 



VNVq 
1 



^{CDC T )- pf {VEC T )^{ yi ). (C4) 



In this way we express antinucleon to pion TDAs (respectively, nucleon DAs) through 
nucleon to pion TDAs (respectively, antinucleon DAs). 

Now we proceed with the generic contributions to the amplitude. Each of three lines 
of the w-channel amplitude graph depicted in Fig. corresponds to a matrix in the Dirac 
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FIG. 10: Typical graphs corresponding to NN — > r y*(q)M(pM)- Left panel: -u-channel 
regime. Right panel: t-channel regime. The crosses represent the virtual photon 

vertex. 

space. We denote these matrices as (Li)p P , (L 2 )f T and (L 3 ) xx . Up to the obvious changes 
due to pn <^ Pni the matrices corresponding to three lines of the t-channel graph are 
{^i)ppi (-^Drr and (Lj ) xx . 

With these notations the spinor structure of the amplitude corresponding to the it- 
channel graph reads 

S u = ^(^)(F(p^)£;) x (L 1 ) p - p (L 2 )^( J L3)xx^pr(Piv, A U )(BU( PN )) X . (C5) 
The spinor structure of the amplitude corresponding to the i-channel graph then reads: 



■If 



Vn 



D^ N ){V(p N )E) x {tfLlC) pP {tfLlC) T ,{& 

(C6) 



Given that for any of 21 relevant graphs presented in Table I of Ref. [13] Li, L 2 and L3 
contain an odd number of the Dirac matrices 7^ (and do not contain 75) we conclude that 
S u turns into i]~ 1 S t with the change ( IClj) . Thus we arrive at the conclusion that the C- 
invariance results in a symmetry between the u- and t-channel factorization mechanisms. 
The amplitudes for NN — > 7*7r within the two mechanisms have the forms, 

1 



M x = C 

■> l S N Sfj Q 4 



V(p„, Sjt)t(\) l5 U(p N , s N )l(C, u) 



-— Vfa, s R )e*(\)A^ l5 U(p N , s N )l'(C, u) 



(C7) 



V(p N , s N )t(X) l5 U(p R , s N )l(e, t) 



—^V( Pn , s N )e*(X)A t Tl5 U(p R , s R )T(e, t) 
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with the same functions X and X' (fT5|) standing both in ( 1C8j) and (1C7j) . Although a 
somewhat unusual choice of the Dirac spinors is employed for S f in flC7j) . this does not 
lead to any particular complications for the calculation of unpolarized cross sections. 
Indeed, ( 1C1I) can be used at the cross section level after the summation over nucleons 
spins is performed. As a result, the forward and the backward peaks of the unpolarized 
cross section possess the identical structure. 
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